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Combined  nonlinear  effects  in  two-photon  absorption  chromophores 

at  high  intensities 

R.  L.  Sutherland*3,  D.  G.  McLean3,  M.  C.  Brant3,  J.  E.  Rogersb,  P.  A.  FIeitzc,  and  A.  M.  Urbase 
3Science  Applications  International  Corporation,  Dayton,  OH,  USA  45431 
bUES,  Inc.,  Dayton,  OH,  USA  45432 

cAir  Force  Research  Laboratory  Wright-Pattcrson  Air  Force  Base,  OH,  USA  45433 

ABSTRACT 

Large  two-photon  and  excited  state  absorption  have  been  reported  in  donor-acceptor-substituted  pi-conjugated  mole¬ 
cules.  We  have  performed  detailed  nonlinear  absorption  and  photophysical  measurements  on  a  system  of  AFX  chromo¬ 
phores  and  calculate  the  nonlinear  transmission  based  on  an  effective  three- level  model.  A  numerical  model  that  in¬ 
cludes  far  wing  linear  absorption  has  been  developed  and  compared  with  an  analytical  three-photon  absorption  model. 
The  models  are  in  accordance  and  yield  excellent  agreement  with  experimental  nonlinear  transmission  data  for  0.G2-M 
AFX  solutions  up  to  laser  intensities  ~  1  GW/cml  Concentration  effects  at  this  intensity  become  increasingly  evident. 
We  have  extended  our  modeling  efforts  to  include  some  new  effects  that  may  be  anticipated  in  this  regime,  such  as 
stimulated  scattering,  molecular  interactions,  and  saturation.  Effects  of  chirped  pulses  and  linewidth  of  the  pump  laser 
on  stimulated  scattering  are  included.  We  report  on  our  experimental  observations  of  various  materials  and  discuss  re¬ 
sults  with  respect  to  our  extended  theoretical  models. 

Keywords:  two-photon  absorption,  excited  state  absorption,  nonlinear  transmission,  nonlinear  scattering 

1.  INTRODUCTION 

Symmetric  and  asymmetric  electron-donor/acceptor-substituted,  Tt-conjugated  systems  are  a  major  class  of  enhanced 
two-photon  absorption  (TPA)  materials.  1,2  The  molecular  TPA  cross  section,  a2i  is  often  characterized  by  nonlinear 
transmittance  (NLT)  experiments,  both  in  the  nanosecond  and  femtosecond  regimes.  However,  the  nanosecond  meas¬ 
urements  in  these  materials  typically  yield  values  of  o2  larger  by  more  than  two  orders  of  magnitude.3,4  Excited  state 
absorption  (ESA)  has  been  postulated  to  play  a  role  in  the  nanosecond  measurements,  and  for  this  reason  the  nonlinear 
parameters  have  been  called  effective  TPA  cross  sections,5,6 

TPA  followed  by  ESA  was  observed  as  early  as  1974,  and  a  value  for  the  product  of  the  ESA  cross  section  and  excited 
state  lifetime  of  the  chromophore  was  estimated  based  on  a  rate  equation  analysis.7  Several  authors  have  posited  the 
effects  of  excited  states  on  nonlinear  absorption  and  refraction,  and  these  effects  have  been  convincingly  demonstrated 
in  degenerate  four-wave  mixing  and  Z-scan  experiments.5'11  Evidence  of  ESA  has  been  demonstrated  in  NLT  and  Z- 
scan  measurements  of  organic  materials12'13,  and  effective  (intensity  dependent)  TPA  coefficients  have  been  employed 
to  estimate  ESA  cross  sections.3,11,14  However,  to  our  knowledge  no  one  has  independently  identified  and  characterized 
these  excited  states,  then  used  this  information  to  theoretically  predict  or  model  nanosecond  NLT  measurements.  The 
nature  of  these  excited  states  is  important.  For  example,  properties  of  the  excited  singlet  state  are  significant  for  fluores¬ 
cence  imaging  and  laser  applications,  whereas  the  triplet  state  plays  an  important  role  in  photopolymerization  and 
photodynamic  therapy.  Moreover,  the  intrinsic  (pulse-width  independent)  value  of  the  TPA  cross  section  a2  as  well  as 
the  significance  of  ESA  in  comparison  to  other  potential  nonlinear  mechanisms,  such  as  stimulated  scattering,  self- 
foeusing/defocusing,  and  possibly  two-step  TPA15,  need  to  be  elucidated  for  accurate  modeling  of  nonlinear  absorption 
in  the  nanosecond  regime.13  Our  goal  is  to  ensure  that  all  of  the  material  parameters  necessary  to  model  the  NLT  are 
independently  characterized.  Eliminating  these  other  phenomena  as  potential  sources  of  NLT  simplifies  the  modeling. 

Recently,  He  et  ah  reported  the  observation  of  stimulated  backscattering  in  a  two-photon  absorption  medium,  where  the 
frequency  of  the  stimulated  wave  was  identical  to  the  incident  laser  frequency  (within  the  resolution  of  their  interfer¬ 
ometer),  and  the  small-signal  gain  was  quadratic  in  the  incident  laser  intensity,1617  They  considered  a  stimulated  thermal 
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Rayleigh  scattering  model  based  on  TPA-enhanced  temperature  and  density  fluctuations,  but  ruled  this  out  due  to  (a)  the 
broad  line  width  of  their  pump  laser,  which  would  severely  reduce  the  gain  of  the  stimulated  wave,  and  (b)  the  fact  that 
the  peak  gain  in  this  theory  is  for  an  anti- Stokes-shifted  wave,  contrary  to  their  experimental  results.  They  concluded 
that  the  stimulated  wave  was  due  to  a  Bragg  grating  formed  by  the  superposition  of  the  incident  laser  beam  and  an  elas¬ 
tically  (Rayleigh)  backscattered  wave,  creating  an  index  grating  via  a  TPA-resonance  enhanced  nonlinear  index  coeffi¬ 
cient  n2 ,  i.e.,  a  Kerr  effect.  The  stimulated  backscattered  wave  is  then  due  to  reflection  of  the  laser  from  this  grating. 

Two-beam  coupling  (TBC)  has  a  rich  history  in  nonlinear  optics.  Early  studies  established  two  conditions  for  one-way 
energy  transfer  between  two  electromagnetic  waves  in  a  Kerr  medium;  1)  the  frequencies  must  be  nondegenerate,  and 
2)  the  nonlinearity  of  the  medium  must  have  a  finite  response  time,1*  The  nonlinear  refractive  index  is  usually  assumed 
to  obey  a  simple  Debye  relaxation  model.  Energy  transfer  under  these  conditions  is  always  from  the  high  frequency 
beam  to  the  low  frequency  beam  for  a  positive  Kerr  nonlinearity,  Silberberg  and  Joseph  showed  that  counter- 
propagating  beams  in  such  a  Kerr  medium  can  exhibit  optical  instabilities  and  self-oscillation,18,19  while  Yeh  gave  an 
exact  solution  for  energy  transfer  between  two  co-propagating  beams  including  the  effect  of  linear  absorption.20  In  each 
of  these  studies,  both  beams  had  non-zero  input  values. 


TBC  also  describes  several  stimulated  scattering  phenomena,  such  as  stimulated  Raman  (SRS),  Brillouin  (SBS),  and 
Rayleigh-wing  scattering  (SRWS).  The  Stokes  (down-shifted  frequency)  wave  in  these  phenomena  arises  internally 
from  scattered  light  and  experiences  gain  at  the  expense  of  the  incident  laser  beam.  The  Stokes  frequency  shift  is  a  nor¬ 
mal  mode  of  the  medium  (e.g.,  vibrational)  for  SRS.  an  acoustic  frequency  for  SBS,  and  the  inverse  reorientation  time 
of  an  anisotropic  molecule  for  SRWS.  SRS,  SBS,  and  SRWS  are  resonant  phenomena,  indicating  that  the  nonlinear 
medium  has  a  finite  response  time.  In  each  of  these  stimulated  scattering  effects  the  small-signal  gain  of  the  Stokes 
wave  is  proportional  to  the  laser  intensity.21 


The  model  proposed  by  He  et  at.  would  appear  to 
violate  the  two  conditions  previously  established 
for  one-way  energy  transfer  by  TBC  in  a  Kerr  me¬ 
dium.  The  frequencies  are  degenerate,  and  the  third 
order  susceptibility  related  to  TPA  is  due  to  elec¬ 
tronic  polarization,  which  has  a  virtually  Instanta¬ 
neous  response  for  nanosecond  laser  pulses.  Al¬ 
though,  under  the  condition  of  a  TPA  resonance, 
this  mechanism  may  indeed  have  a  finite  response 
time,  the  intensity  dependence  of  the  Kerr  effect  in 
TBC  is  inconsistent  with  the  experimental  results 
reported  by  He  et  ai. 

We  have  performed  detailed  nonlinear  absorption 
and  photophysical  measurements  on  a  system  of 
AFX  chromophores  and  calculate  the  nonlinear 
transmission  based  on  an  effective  three-level 
model,  A  numerical  model  that  includes  far  wing 
linear  absorption  has  been  developed  and  compared 
with  an  analytical  three-photon  absorption  model. 
The  models  are  in  accordance  and  yield  excellent 
agreement  with  experimental  nonlinear  transmis¬ 
sion  data  for  0,02-M  AFX  solutions  up  to  laser  intensities  -  1  GW/cm2.  Concentration  effects  at  this  intensity  become 
increasingly  evident.  We  have  extended  our  modeling  efforts  to  include  some  new  effects  that  may  be  anticipated  in  this 
regime,  such  as  stimulated  scattering,  molecular  interactions,  and  saturation.  Effects  of  chirped  pulses  and  linewidth  of 
the  pump  laser  on  stimulated  scattering  are  included.  We  report  on  our  experimental  observations  of  various  materials 
and  discuss  results  with  respect  to  our  extended  theoretical  models. 
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2.  MATERIALS 


The  chromophores  used  in  this  study  are  part  of  a  class 
of  TPA  materials  (designated  AFX)  having  a  design 
based  on  a  push-pull  charge-transfer  model  and  multi¬ 
dimensional  conjugation  motif6  Examples  of  these 
molecules  are  shown  in  Fig.  1 .  AF24G  is  a  linear  D-tt- 
A  chromophore,  where  the  donor  (D)  is  diphenyl- 
amino  and  the  acceptor  (A)  is  benzothiazole.6  AF350 
is  a  three-arm  octupolar  molecule  (D-A3)  with  a  single 
triarylamino  group  serving  as  the  electron-rich  hub, 
dialkylfluorenyl  bridges,  and  three  7r-electron  deficient 
benzothiazoles.422  Except  for  one  less  phenyl  ring, 
AF24Q  is  essentially  one  arm  of  AF350.  AF455  is  an 
A-D3  type  octupolar  molecule  with  1,3,5-triazine  as 
the  hub  moiety  and  dipheny famines  as  the  Ti-electron 
donating  end-groups,4  These  molecules  were  synthe¬ 
sized  and  purified  in  the  Polymer  Branch  of  the  Air 
Force  Materials  Directorate  (AFRL/MLBP)  as  de¬ 
scribed  in  References  4,  6,  and  22*  Solutions  of  each 
were  prepared  in  tetrahydrofuran  (TFIF)  for  both  pho¬ 
tophysical  characterization  and  NLT  experiments* 

A  typical  linear  absorption  spectrum  is  shown  in  Fig* 
2,  The  Xmax  is  in  the  UV,  with  absorption  band  edges 
ranging  from  ~  450  nm  to  ~  470  nm,  In  each  case 
there  is  very  little  absorption  (o0  <  10'2]  cm2,  see  Fig. 
2b)  in  the  2*5-1  *4  eV  (500-900  nm)  region*  The  energy 
difference  between  the  absorption  band  edge  and  a 
laser  line  at  800  nm  is  -  1.5  eY.  The  purity  should  be 
at  least  98%  for  glassy  materials  (AF-455)  and  99% 
pure  for  crystalline  compounds  (AF-240,  350)  based 
on  elemental  analysis  together  with  spectroscopic 
methods  and  melting  point  determination.  There  is 
only  a  single  peak  detected  by  liquid  chromatography 
for  both  AF-350  and  AF-455.  An  impurity  with  a  peak 
in  the  500-900  nm  range  with  a  peak  extinction  coeffi¬ 
cient  of  >10,000  M^cm"1  would  be  observed  at  a  frac¬ 
tional  concentration  of  lOppm  or  less.  We  conclude 
that  it.  is  unlikely  that  there  is  significant  impurity  ab¬ 
sorption  in  the  samples. 

3.  THEORETICAL  MODELS 


Photon  Energy  (eV) 


3.1.  Effective  TPA  and  3PA  Models 


Fig,  2,  Linear  absorption  spectra  (ground  state  cross  section  o0) 
of  AF240,  AF350,  and  AF455  solutions  in  THF(eides),  The  solid 
lines  are  a  Voigt  function  fit  to  the  data. 


We  seek  a  simple  analytical  model  that  can  adequately 
explain  and  reliably  predict,  within  a  reasonable  de¬ 
gree  of  approximation,  the  nanosecond  NLT  based  on 

measurable  properties  of  the  chromophores  and  the  laser  pulse.  This  model  is  similar  to  that  used  for  reverse  saturable 
absorption  media,  but  with  TPA  from  the  ground  state.23  We  consider  ESA  from  the  lowest  lying  singlet  and  triplet 
states.  Two-photon- induced  ESA  is  a  three-photon  process,  although  three  photons  are  not  absorbed  simultaneously  due 
to  the  finite  lifetimes  of  the  excited  states.  It  is  of  interest,  nevertheless,  to  see  how  useful  a  three-photon  absorption 
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(3 PA)  model  is  in  explaining  the  data.  In  the  following  we  present  the  conditions  and  approximations  for  employing 
such  a  model. 

With  the  above  assumptions,  we  have  for  the  intensity  /, 

~  =  +  (I) 

oz 

where 


T  = 


=TlaS1  +fe-Tl)pTaTi. 


(2) 

(3) 


and  p  =  g2jV  is  the  TP  A  coefficient.  We  have  thus  modeled  the  system  as  a  two-photon  resonant,  three-photon  absorp¬ 
tion  medium,  with  an  effective  intermediate  state  that  represents  a  time-averaged  combination  of  the  and  T]  states. 
Several  authors  have  taken  the  quantity  in  parentheses  in  Eq.  (1)  to  represent  an  effective  TP  A  coefficient  [3^,  from 
which  an  effective  TPA  cross  section  can  be  derived.  We  will  consider  both  effective  TPA  and  3PA  models  in  analyz¬ 
ing  the  experimental  data.  The  quantity  r\  is  the  pulse-averaged,  relative  weight  of  the  Si  state  contribution  to  the  effec¬ 
tive  ESA  cross  section.  When  xp  «  ts,  r\  &  lA  and  ESA  is  essentially  from  the  state  only.  On  the  other  hand,  when  xp 
»  Ts,  *  0  and  all  ESA  is  effectively  from  the  T|  state.  When  xp  «  xs,  the  growth  of  the  St  population  is  linear,  and 
the  average  growth  time  over  a  rectangular  pulse  is  just  ip/2,  The  maximum  relative  weight  of  VA  for  either  state  reflects 
the  fact  that  in  either  time  regime  the  growth  rate  of  the  respective  excited  state  population  (singlet  or  triplet)  is  ap¬ 
proximately  constant.  The  average  population  density  over  the  entire  pulse  is  then  just  one-half  of  the  maximum  density 
obtained  at  the  end  of  the  pulse,  iVsl  max  =  (pT  NJ}  max  =  tpG2Ar/peak/2/k][) .  In  an  intermediate  regime  where  xT  »  xp 
and  xp  -  xs,  tj  <  lA  means  that  the  Sj  population  at  the  end  of  the  pulse  is  smaller  than  Nsl  because  of  con¬ 

tinuous  transitions  to  the  So  and  Ts  states  as  Si  is  being  pumped,  analogous  to  a  leaky  capacitor. 

Equation  (1)  can  be  solved  analytically,  but  the  resulting  expression  is  transcendental.  However,  if  we  take  (3  «  y /,  the 
result  reduces  to  a  simple  closed  form  identical  to  that  for  three-photon  absorption.  Assuming  now  Gaussian  spatial  and 
temporal  profiles  for  the  incident  intensity,  we  integrate  the  resulting  expression  over  space  and  time  to  obtain  the  en¬ 
ergy  transm  ittance  for  a  sample  of  thickness  d\ 


T  = 


JxPo  }  ' 
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V1  +  pi  exp(-2jf2J  +  p0  exp(- 


*2) 
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where  To  is  the  net  linear  transmittance  from  air  into  the  medium  due  to  Fresnel  losses  only  (including  all  dielectric  in¬ 
terfaces),  po  =  (2yro2/oV)l/2,  and  /0  is  the  incident  peak,  on-axis  intensity.24  The  parameters  determining  y  can  be  found 
from  femtosecond,  picosecond,  and  nanosecond  photophysical  measurements  and  standard  TPA  measurements  in  the 
femtosecond  regime. 


3,2*  Numerical  ESA  Model 


The  numerical  model  calculates  the  spatial,  time,  and  radial  dependence  of  three  state  populations;  all  of  the  transitions 
between  these  states;  the  bimolecular  processes  of  trip  let- triplet  annihilation,  ground  state  quenching,  and  oxygen 
quenching.  The  radiation  transport  equation  is  then, 

31  2 

~  =  ”°SoAV  “  °2^$0^  '  “  GSI^S1^  —  aTt^Tl^ 

&  ■  (5) 


The  population  rate  equations  are, 


4 


dN. 


so 


dt  h(fi 


f  o3f) 

°so  K^so  + 

1 

— i 

_ 

( 

^si  + 

l  2  J 

k  Ts  J 

k 

K. 

”  •"  _  J^  TI  +  KST^SO  +  KOT^O 

r  t  Z 


AT. 


Tl 


5M 


SI  _  ^  _  ,  G2^ 


3/  ftco 


= —  I  aso  +  — |iV™  - 


N. 


si 


"so 


■jV 


Tl 


dNTl  _  <PX  „ 
dt  t 


SI 


+  NTl  +  K  ST  N$q  +  k0T7Vq 


V  rT  2 


TV- 


Tl 


3A( 


dt 


—  KOT  ^O^TI  +kotYo  JVo) 


NS0=N-NS]-Nll 


(6a) 

(6b) 

(6c) 

m 

(6e) 


where  Krr,  KsTi  KqT,  TV0,  ko  are  the  triplet-triplet  annihilation  rate,  the  ground  slate  self  quenching  rate  of  the  triplet 
state,  the  oxygen  quenching  rate  of  the  triplet  state,  the  oxygen  concentration,  and  the  oxygen  singlet  state  relaxation 
rate  respectively.  The  rate  equations  are  solved  using  a  stiff  Runge-Kutta  algorithm  supplied  in  MathCad.  This  requires 
the  Jacobian  to  be  supplied  which  was  derived.  This  is  solved  at  every  time  step.  These  populations  are  used  to  calculate 
the  change  in  the  intensity  which  is  then  rescaled  assuming  Gaussian  beam  propagation  at  each  z  step.  A  range  of  input 
intensities  are  calculated  and  used  to  calculate  both  the  NLT  variation  and  the  radial  dependence.  A  Gaussian  time  pro¬ 
file  is  used  since  the  measured  laser  profile  nearly  matches  a  Gaussian,  The  radial  and  temporal  integrations  are  done  to 
give  transmitted  pulse  energy  and  compared  to  the  measured  data.  A  normal  reflection  Fresnel  loss  is  imposed  on  the 
incident  intensity  and  again  on  the  output  intensity.  A  refractive  index  of  1.4535  is  used  for  fused  silica  at  800  nm  and 
1.51 1  for  BK7  at  800  nm.  The  Fresnel  reflection  at  the  glass  liquid  interface  is  found  to  be  negligible  and  is  not  in¬ 
cluded. 


3,3,  Stimulated  Scattering  Model 

Consider  a  field  with  a  time  dependent  amplitude  and  phase,  ^(/)exp[- />(*}],  A  Taylor  series  expansion  of  the 

phase  yields  (ignoring  a  constant  term)  #>(/)=  cot  +  bt2  +  ,  where  w  is  the  central  frequency  of  the  wave,  and  b  is  a  lin¬ 

ear  chirp  coefficient.  For  simplicity,  I  will  ignore  higher  order  terms.  Co-propagating  TBC  with  chirped  pulses  has  been 
examined  in  Kerr  media  with  a  finite  response  timer5-  However,  none  of  these  previous  studies  have  considered  the 
interaction  of  counter-propagating  waves  through  TPA-populated  excited  states. 

Let  the  total  field  in  a  medium  of  length  d  be  described  by 

E(z ,  t )  =  A L (z,  t  +  r) exp j  i[kz  -cot  - b(t  +  r)2  ]  } 

+  {z,  t  -  r)exp  j  -cot-b{t  -  r)2  ] }+  c,c, 

where  k  =  ncdc,  n  is  the  linear  refractive  index,  and  2r-  2 n(d -  z)!c  is  the  relative  time  delay  at  position  z  and  time  i 
between  the  forward  propagating  laser  wave  (L)  and  the  backward  propagating  scattered  wave  (S).  I  assume  that  the 
scattered  wave  originates  at  z-d  by  some  elastic  scattering  process,  so  As(dj)- ^rjAL(d,t)  where  tj  is  a  constant  «  L 
Consequently,  the  scattered  wave  has  the  same  spectral  composition  as  the  incident  laser  wave.  However,  for  z  ^  d  a 
lower  frequency  part  of  the  scattered  wave  is  always  interacting  w  ith  a  higher  frequency  part  of  the  incident  wave.  The 
total  polarization  of  the  medium  is  gi  ven  by 

P-e  oUl>  +  +3/jt(3>(£2})e  (S) 

where  is  the  ^-th  order  susceptibility,  and  the  angular  brackets  indicate  an  average  over  a  time  longer  than  an  optical 
period  but  shorter  than  (2br)[ ,  We  have  assumed  that  TP  A  produces  a  single  excited  state  of  number  density  Ne  «  N , 
the  total  number  density  of  the  nonlinear  ehromophore.  We  have  also  assumed  an  isotropic  medium  w  ith  linearly  polar- 
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ized  light  for  simplicity,  and  both  Ay{1)  =  %{P  -  and  are  complex  quantities  (g  signifies  the  ground  state  of  the 
medium).  The  excited  state  decays  back  to  the  ground  state  with  a  time  constant  Tey  so  Ne  obeys  the  following  kinetic 
equation: 

tr2A72  at, 

dt  2  tm  Te  K  } 

where  cr,  =  p!N  is  the  TP  A  cross  section  (/?  is  the  TP  A  coefficient),  and  /(z,/)  =  2e0nc^E2(z,tfj  is  the  total  intensity.  Let 
the  amplitudes  be  slowly  varying  in  time  compared  to  Te*  Equation  (9)  can  then  be  integrated  to  yield 

NJN  =  C0  +  {C,  exp|/(2£z  -  4 bzt)]+  C2  exp[i{4kz  -  8 h  * )]+  cjc) ,  ( 1  Oa) 


c0=^h+/j+444), 

2  hm 

(10b) 

2 hm  \  1  -  iAb  tTc  J 

(10c) 

c  (2 e0ncA,  As)2 

2  2 hm  1  -  iSb  rTe  ’ 

(10d) 

where  £q  is  the  free-space  permittivity.  We  will  make  the  assumption  that  (4 bzTe)2  «  1. 

Following  a  procedure  directly  analogous  to  that  of  Yeh:o  and  Boyd:i  for  non-degenerate  TBC,  Eqs.  (10a)-(10d)  are 
substituted  into  Eq,  (8),  and  then  both  Eqs.  (7)  and  (8)  are  substituted  into  Maxwell's  wave  equation  in  the  slowly  vary¬ 
ing  amplitude  approximation.  Matching  up  synchronous  terms,  the  following  coupled-wave  equations  for  the  laser  and 
back  scattered  intensities  are  derived: 

=  -*( l  -  4^4  +  4)44  -  r.jM  +  3/i  +  644 K  -  p(!i  +  24)4  Ola) 

~  =  -g(!  -  444  +  4  )44  +  r+tyl  +  /!  +  644  K  +  P(2It .  +  4)4  (lib) 


where  g  and  y^are  the  backward  wave  gain  and  effective  three- photon  absorption  (3 PA)  coefficients,  respectively,  with 


%bTem\xfd 


(12) 


~ 


N&a 


(13) 


where  A^11  =  A  a  =  ae  -a  is  the  difference  between  the  linear  absorption  cross  sections  of  the  excited  and 

ground  states,  and  I mt  =  {lhcolcr1T^u  is  the  two-photon  saturation  intensity.  Note  that  g  =  0  when  b  =  0  (no  chirp). 
Hence,  without  chirp  there  can  be  no  growth  of  the  backward  scattered  wave.  The  gain  g  also  carries  the  sign  of  b  (posi¬ 
tive  or  negative  chirp).  For  a  negatively  chirped  pulse,  the  scattered  wave  will  be  attenuated.  For  the  rest  of  the  paper, 
we  will  assume  that  b  >  0. 


4*  EXPERIMENTAL 

Micromolar  solutions  were  prepared  for  conventional  photophysical  measurements  (lifetimes,  excited  state  cross  sec¬ 
tions,  and  quantum  yields).  Details  of  these  experiments  can  be  found  elsewhere."3  All  solution  samples  for  nanosecond 
NLT  measurements  had  concentrations  of  0,02  M  (mol/L)  and  were  placed  in  1-mm  glass  or  fused  silica  cuvettes.  In¬ 
trinsic  g2  values  were  obtained  from  independent  femtosecond  measurements.23  25 
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Nanosecond  nonlinear  transmittance  measurements  were  performed  with  a  Nd:YAG  pumped  optical  parametric  oscilla¬ 
tor  (OPO)  tuned  from  660  to  880  nm.  The  pulse  was  Gaussian  shaped  with  =  3.2  ns.  The  beam  was  focused  with  an  f 
50-cm  lens  into  the  sample.  Over  the  length  of  the  sample  {1  mm)  the  beam  was  essentially  collimated.  The  beam 
shape  was  slightly  elliptical,  with  a  geometric-mean  1/e2  radius  w  »  18,4  -  18.9  pm,  assuming  an  approximately  Gaus¬ 
sian  beam  shape.  7 he  energy  was  varied,  and  incident  and  transmitted  energies  vvere  measured  with  energy  meters.  To 
rule  out  the  effects  of  self-focusing/defocusing,  a  large-area  (~  1  cm:)  detector  was  placed  near  the  exit  of  the  sample  to 
collect  all  transmitted  energy.  We  also  looked  for  stimulated  backscattering  by  rotating  the  sample  slightly  to  avoid 
Fresnel  reflected  light  and  measuring  1 80°-scattered  light  with  an  energy  meter. 

Intermolecular  quenching  properties  were  examined.  The  ground  state  concentration  was  varied  from  10  pM  to  20  mM, 
and  the  triplet  state  lifetime  measured.  At  -  10  pM  the  lifetime  of  AF350  was  182  ps  under  deoxygenated  conditions, 
and  the  lifetime  was  96ps  at  20  mM.  A  bimolecular  self-quenching  rate  constant  of  2  x  I02  M  's"1  was  measured  The 
NLA  experiments  were  performed  under  air-saturated  conditions,  and  this  small  amount  of  quenching  will  make  no 
difference  since  it  is  in  competition  with  oxygen  quenching.  Similarly,  AF350  run  at  low  and  high  energy  excitation 
conditions  produced  no  measurable  change  in  the  decay  rate  attributable  to  triplet-triplet  annihilation.  This  is  not  unex¬ 
pected  due  to  the  low  triplet  yield  of  these  materials.  Due  to  the  similarity  of  these  materials  in  structure  and  known 
triplet  yield,  we  assume  that  these  intermolecular  processes  do  not  play  a  significant  role  at  this  concentration. 

5.  RESULTS  AND  DISCUSSION 

A  common  method  in  the  literature  for  screening  TPA  materials  in  the  nanosecond  regime  is  to  fit  the  NLT  data  to  an 
effective  TPA  transmittance  (see  Eq.  (1))  and  extract  an  effective  TPA  cross  section.  We  show  the  results  of  such  fits 
for  two  series  of  dipolar  (AF240  and  AF  270),  quadrupolar  {AF287  and  AF295),  and  octupolar  (AF380  and  AF350) 
molecules  in  Table  1.  Another  series  of  chromophores  have  been  measured  at  a  variety  of  wavelengths,  and  their  effec¬ 
tive  TPA  cross  sections  are  given  in  Table  2.  Such  data  are  useful  for  observing  trends,  in  molecular  structure  and  spec¬ 
trally,  but  do  not  yield  sufficient  information  to  give  directions  for  further  improvement.  For  example,  a  large  effective 
o2  could  be  due  to  a  large  intrinsic  o2,  a  large  singlet  and/or  triplet  cross  section,  and/or  a  large  triplet  yield.  Flence, 
these  NLT  measurements  must  be  supplemented  with  photophysical  measurements.  Both  effective  3PA  and  numerical 
ESA  models  may  then  be  applied  to  ascertain  the  quality  of  the  model  in  reproducing  the  NLT  results. 

Table  1.  Effective  TPA  Cross  Sections  of  a  Series  of  AFX  Chromophores 


Material 

File 

Effa2 

( 1  O'20  cm4/GW) 

AF24Q 

NLA  130 

50 

AF287 

NLA  134 

99 

AF3S0 

NLA  131 

114 

AF270 

NLA  132 

29 

AF295 

NLA  136 

73 

AF350 

NLA  137 

139 

7 


Table  2.  Spectral  Effective  TPA  Cross  Sections  for  Four  Chromophores 


Chromophore 

File 

EfF  o2 

{10’2°  cm4/GW) 

X  (nm) 

AF445 

NLA0841 

44 

660 

NLA0771 

103 

740 

N  LAO  8  04 

29 

800 

N  LAOS  16 

8 

880 

El-BTF 

N  LAO  8  42 

195 

660 

NLA0756 

145 

740 

NLA0801 

80 

800 

N  LAOS  17 

0 

880 

AF380-1  18 

NLA0847 

67 

660 

NLA0761 

117 

740 

NLA0807 

69 

800 

N  LAO  8  22 

100 

880 

IR-2 

NLA0850 

114 

660 

NLA0776 

98 

740 

NLA0809 

103 

800 

NLA0825 

185 

880 

Results  of  NLT  measurements  for  0.02-M  solutions 
of  AF455,  AF350,  and  AF240  in  THF  are  given  in 
Figs,  3-5,  respectively.  Using  sets  of  intrinsic  cr2 
values,  we  calculated  the  NLT  for  each  material  and 
compare  these  calculations  to  the  data  in  Figs.  3-5. 

Only  experimental  data  were  used  in  the  calcula¬ 
tions;  there  were  no  adjustable  parameters. 

In  Fig.  3  there  is  good  agreement  between  theory  and 
experiment  when  the  value  of  cr2  -0.51  x  10'20 
cm4/GW  is  used.  This  value  was  obtained  in  a  fem¬ 
tosecond  NLT  experiment30  for  a  0.02-M  solution  of 
AF455  in  THF,  identical  to  the  sample  studied  here. 

Although  the  wavelength  used  in  that  experiment 
was  790  nm,  the  TPA  spectrum4  indicates  that  g2 
does  not  differ  significantly  at  800  nm,  considering 
the  experimental  uncertainty  of  ±15%. 

Figure  3  also  compares  the  transmittance  in  AF455 
calculated  for  two-photon  induced  ESA  (effective 
3 PA  analytical  model)  with  that  due  to  TPA  alone 
(numerical),  and  reverse  saturable  absorption  (RSA) 
which  includes  the  ground  state  absorption  and  ex¬ 
cited  state  absorption.  Obviously,  ESA  is  the  domi¬ 
nant  loss  mechanism.  The  inclusion  of  ground  state 
absorption  gives  rise  to  significant  nonlinear  absorp¬ 
tion  and  fits  the  data  more  closely  in  the  region  near  <  10  pJ  where  these  models  are  the  most  reliable. 


Fig,  3,  Experimental  (circles)  and  theoretical  (curves)  nanosecond 
nonlinear  transmittance  as  a  function  oflaser  pulse  energy  for  a  0  02- 
M  THF  solution  of  AF455  at  800  nm.  The  curves  are  for  the  analyti¬ 
cal  effective  three-photon  absorption  model,  the  numerical  model 
with  only  TPA  driving  excitation,  with  only  the  ground  state  driving 
excitation  (labeled  RSA)  and  for  the  complete  system. 
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Fig.  4.  The  curves  are  as  described  in  Figure  1 0  but  for  a  Q.Q2-M  THF 
solution  of  AF350  at  820  nm 


We  note  for  all  materials  that  there  is  some  departure  of  Fi8-  5-  The  curves  are  as  described  n  Figure  10  but  for  a  0.02-M  THF 
experiment  from  theory  in  the  low  energy  regime  (~  10  solution  of  AF240  at  820  nm,. 

pjf  It  is  interesting  that  an  effective  TPA  theory  fits  the  data  better  in  this  regime,  indicating  that  some  mechanism  that 
has  a  quadratic  intensity  dependence  may  be  at  play.  He  et  al,37  have  observed  this  kind  of  dependence  and  have  veri¬ 
fied  in  their  case  that  this  is  due  to  a  stimulated  scattering  phenomenon.  We  consider  this  next. 


Let  us  examine  Eqs.  (1  la)  and  (I  lb)  in  the  case  when  Is  «  h  55  constant.  It  can  be  seen  that  the  backward  stimulated 
wave  will  experience  exponential  growth  when  IL  >  This  defines  the  threshold  condition  for  stimulated 

backscattering  when  there  is  no  linear  absorption.  From  Eq.  (lib),  the  small  signal  gain  G ,  given  in  terms  of 
Us=Is(o)-Is(d)>  is  G  =  AIS  fls  (</) oc  if  for  input  intensities  where  the  gain  dominates  TPA.  Also,  the  reflectance  of 
the  scattered  wave  is  defined  by  R  =  IS{ o)/!L  (o),  and  the  change  M  =  R-rj  =  &Is/U  * h(d)h  ^ 0!s{d)j} L  .  These  re¬ 
sults  are  in  agreement  with  the  data  presented  by  He  et  al.  for  the  measured  small-signal  gain  and  reflectance  in  a  0.01- 
M  solution  of  PRL  802  in  THF.t7 


Under  conditions  where  the  nonlinear  absorption  terms  (0  and  yeJ$  can  be  ignored,  Eqs.  (1  la)  and  (11b)  can  be  solved 
analytically.  The  result  can  be  expressed  as 


R 

[{l-/fXl-JR  +  377o)+2?7o2] 

1/2 

fl-rt  +  2/?^ 

(l  +  r)2  exp(f) 

v  l-«  +  7o  j 

(14) 


where  rj0  -  and 


T  =  \{\-Rfgll{v>)d. 


(15) 


Note  that  in  general  rj0  *  ?h  although  when  the  gain  is  sufficiently  small  rj o  ^  77.  Also,  when  the  gain  is  not  too  large  so 
that  (1  -  R)»  tj o,  Eq.  (14)  can  be  simplified  to  the  following  approximation: 


/;(i  -  /?) 

(i  +  /?)2  exp(r) 


06) 


It  is  interesting  to  compare  this  result  to  the  case  of  SBS,  for  which  the  denominator  of  Eq,  (16)  becomes  exp(r)-fl 
with  F  -*■  (l-^)gB/,  (o)d  ,  where  gB  is  the  Brillouin  gain  coefficient.21  Thus,  in  the  case  of  negligible  loss  by  absorption, 
the  backscattered  reflectance  as  a  function  of  incident  laser  intensity  will  look  similar  to  the  Brillouin  reflectance,  A  plot 
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of  the  reflectance  given  by  Eq.  (16)  for  a  constant  input  0),  over  a  range  off  where  the  approximation  is  valid, 

is  shown  in  Fig.  6. 

Numerical  solutions  of  Eqs.  ( 11a)  and  (1 1b)  are  also  given  in  Fig.  6+  These  illustrate  the  agreement  of  the  approximation 
given  by  Eq.  (16)  with  the  exact  result,  and  the  effect  of  including  the  nonlinear  absorption  terms.  Figure  7  gives  /s{0)  as 
a  function  of  4(0).  Parameters  have  been  adjusted  to  yield  results  close  to  the  experimental  data  of  He  et  al,36’17  The 
coefficients  used  for  the  solid  curve  in  Fig.  7  are  /?  =  9,46  cm/GW  (the  value  quoted  in  Ref  17),  rj  =  0.04,  g  =  1800 
cm3 /GW2,  and  120  cmVoW2,  with  d  =  1  cm.  In  Fig.  8  the  corresponding  plot  for  the  laser  transmittance  is  given. 
Here  the  transmittance  is  the  ratio  of  output  power  to  input  power.  To  compute  this,  we  assumed  a  Gaussian  dependence 
for  the  input  intensity  and  then  integrated  the  output  intensity  over  the  area  of  the  cylindrical ly  symmetric  beam.  The 
result  is  compared  with  what  would  be  expected  for  pure  TP  A  (no  TBC  or  ESA).  The  departure  from  pure  TP  A  in  this 
case  is  due  primarily  to  energy  transfer  from  the  laser  to  the  backscattered  beam.  Figure  7  also  shows  a  somewhat  dif¬ 
ferent  result  (dashed  curve)  yielding  comparable  backscattered  intensity.  Often  the  effective  TP  A  cross  section  meas¬ 
ured  in  the  nanosecond  regime  is  a  factor  -  100  larger  than  the  intrinsic  cross  section  {usually  measured  in  the  femto¬ 
second  regime),31  Thus,  for  these  calculations  we  chose  ft  =  0.09  cm/GW.  The  other  parameters  are  t]  =  0.02,  g  =  2860 
cm3/GW2,  and  yeff=  300  cm3/GW\  TP  A  is  low  in  this  case,  but  the  loss  due  to  ESA  significantly  slows  down  the  growth 
of  the  backscattered  wave  at  higher  intensities. 

To  get  an  order  of  magnitude  of  the  numbers  involved,  consider  the  experiment  of  He  et  al. 16,17  and  the  results  shown  in 
Figs.  7  and  8  for  /?  =  9.46  cm/GW.  For  a  laser  wavelength  of  532  nm  and  Te  ~  1  ns,531  Isat  ~  700  MW/cm2.  For  a  se¬ 
verely  chirped  pulse  the  laser  linewidth  is  A ^  -  2 brL  where  zL  is  the  laser  pulse  width.  For  24  GHz  (0,8  crrf 

l),  Ti~  10  ns,  and  N  =  6  x  10ls  cm'3  (0.01 -M  concentration),  Eqs,  (12)  and  (13)  yield  A^1*  -  4  x  lO’3  and  Act-  1  x  IQ'17 
cm2.  We  note  that  this  numerical  example  contravenes  the  approximation  made  earlier  for  which  (4 brTe)2  is  small  com¬ 
pared  to  1.  To  account  for  deviations  due  to  this,  the  gain  and  ESA  terms  in  Eqs.  (11a)  and  (lib)  would  need  to  be 
modified  by  the  inclusion  of  the  factor  [1  +  (46rre)2]'3.  For  example,  in  the  negligible  nonlinear  absorption  case,  Eq. 
(16)  would  be  multiplied  by  a  factor  ln(l  +  a)/ a,  where  a  =  {AbTend!cfy  which  is  - 1  for  a  «  1.  In  the  present  numeri¬ 
cal  example,  this  factor  is  -  0.5.  Consequently,  the  estimates  for  Aj^l}  and  Act  should  be  increased  by  a  factor  -  2, 
These  values  yield  approximate  agreement  with  photophysical  measurements,  and  thus  suggests  that  this  stimulated 
scattering  with  a  quadratic  intensity  dependence  may  be  contributing  to  the  overall  shape  of  the  NLT  curve. 


Fig.  6.  Backscattered  reflectance  as  a  function  of  exponential  gain  factor  for  various  values  of  nonlinear  absorption  coefficients. 
I^dyiL(Q)  =  10 2  in  all  cases.  Circles  give  values  for  R  calculated  by  the  approximation  in  Eq.  (10), 
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Fig  7  Backscattered  intensity  as  a  function  of  incident  laser  intensity. 


Fig.  8.  Nonlinear  transmittance  of  the  incident  laser  as  a  function  of  incident  intensity. 

Notice  that  Re(^(3)),  or  n2 ,  does  not  appear  in  the  gain  coefficient  of  Eq.  (12).  The  reason  for  this  is  that  was  as¬ 
sumed  to  have  an  instantaneous  response  [  Im(^(3))oc  /?  ].  There  is  no  two-beam  coupling  in  a  Kerr  medium  with  an  in¬ 
stantaneous  response.2021  It  is  possible  however,  under  TEA  resonant  conditions,  that  n2  could  have  a  finite  response 
time  (he.,  ^3)  is  complex  with  a  finite  damping  coefficient21).  However,  if  this  is  included  in  the  development  of  the 
coupled- wave  intensity  equations,  it  would  lead  to  a  term  in  the  gain  G  that  is  proportional  to  lLf  not  /l ,  which  would 
be  inconsistent  with  the  experimental  results  of  He  et  al.1"  It  is  also  quite  likely  that  the  Kerr  refractive  term  would  be 
much  smaller  than  the  term  proportional  to  . 

The  central  feature  of  this  theory  is  contained  in  the  term  4 brTe  =  &&T€y  where  Am  -  4br  represents  the  instantaneous 
difference  in  the  frequency  between  the  forward  and  backward  propagating  waves  at  position  z  in  the  medium.  Am  var¬ 
ies  continuously  from  0  at  z  =  d  to  a  maximum  of  4 hndic  at  z  =  0,  The  interference  of  the  forward  and  backward  propa¬ 
gating  waves  sets  up  an  interference  pattern  that  is  traveling  to  the  left  if  b  >  0,  This  intensity  pattern  forms  a  population 
grating  via  TPA,  which  lags  behind  the  interference  pattern  due  to  the  finite  response  time  Te.  By  examining  Eqs.  (10c), 
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(10a),  and  (8),  we  see  that  Att>re  *  0  results  in  a  contribution  to  the  imaginary  part  of  the  total  refractive  index  of  the 
medium  due  to  the  presence  of  the  two  waves.  As  Boyd  has  observed  in  the  context  of  nondegenerate  TBC  in  a  Kerr 
medium,  the  only  way  to  obtain  a  complex  refractive  index  {apart  from  pure  absorption),  and  thus  to  achieve  energy 
coupling  between  the  two  waves,  is  that  the  product  Aoj7^  not  vanish,21  In  the  present  theory,  this  requires  a  non-zero 
chirp  ( b  *  0).  This  can  also  be  explained  by  considering  the  phase  difference  A<p  =  4brTe  between  the  polarization  and 
the  field.  The  rate  of  energy  transfer  per  unit  volume  to  or  from  an  electric  field  £  by  a  polarization  P  is 
2Re[/^£[[£|exp(/A^)],  which  is  0  if  A (p=  0.32  Thus,  although  backscattering  by  an  unchirped  wave  could  also  lead  to  a 
population  grating,  the  grating  would  not  transfer  energy  to  the  scattered  wave  because  the  field  and  polarization  would 
be  in  phase.  There  would  thus  be  no  increase  in  the  index  modulation  as  the  incident  intensity  increases  and  no  expo¬ 
nential  growth  of  the  scattered  wave,  i,e.,  no  stimulated  scattering.  This  has  an  analogy  in  SBS  and  SRS.  In  both  cases, 
the  gain  is  related  to  the  imaginary  part  of  the  nonlinear  susceptibility,  which  vanishes  when  the  difference  between  the 
laser  and  scattered  wave  frequencies  shrinks  to  zero.  In  addition,  the  scattered  wave  in  the  present  theory  is  attenuated 
for  a  negative  chirp,  analogous  to  the  attenuation  of  the  anti-Stokes  wave  in  SBS  and  SRS.21 

It  should  be  noted  at  this  stage  that  linear  chirp  is  not  a  unique  requirement  for  the  energy  transfer  between  incident  la¬ 
ser  and  elastically  backscattered  waves.  Higher  order  chirp  terms  in  the  nonlinear  time  dependent  phase  have  been  ig¬ 
nored  for  simplicity  but  will  make  additional  contributions  to  the  beam  coupling.  Another  potential  mechanism  involves 
multimode  beams.  In  this  case,  though,  the  spectral  nature  of  the  backscattered  beam  will  not  match  that  of  the  incident 
beam.  Take  the  simple  case  where  the  incident  laser  wave  consists  of  two  modes:  uh  and  m +  Sa>>  (So  «  o*>). 
Employing  the  mechanism  involving  a  TPA-populated  excited  state  described  above,  there  will  be  energy  transfer  from 
the  6> i  mode  of  the  incident  laser  beam  to  the  mode  of  the  backscattered  beam,21  Likewise,  the  a>\  mode  of  the  back- 
scattered  wave  will  yield  its  energy  to  the  mode  of  the  laser  wave.  The  backscattered  wave  will  thus  be  single-mode 
and  its  spectrum  conspicuously  different  from  the  incident  laser.  However,  if  the  modes  in  the  incident  laser  beam  are 
equally  chirped,  both  modes  of  the  backscattered  wave  wil!  be  amplified,  but  not  equally.  When  the  laser  linewidth  is 
determined  primarily  by  the  linear  chirp  (2 brL  »  So?),  the  spectrum  of  the  backscattered  wave  will  superficially  resem¬ 
ble  that  of  the  incident  laser,  but  the  energy  distribution  amongst  the  modes  will  differ.  This  will  generally  be  the  case 
also  when  the  number  of  modes  is  greater  than  two.  Scattering  due  to  multimode  effects  have  not  been  studied  to  any 
large  extent.  We  plan  to  further  examine  this  both  theoretically  and  experimentally. 

6.  CONCLUSIONS 

In  summary,  we  have  measured  the  excited  state  properties  of  donor-acceptor  push-pull  charge-transfer  chromophores 
and  modeled  the  nanosecond  NLT  as  an  effective  three-photon  absorption  process  combined  with  stimulated  scattering 
at  high  intensities.  All  model  parameters,  including  intrinsic  TPA  cross  sections,  have  been  measured  independently, 
and  the  model  agrees  very  well  with  much  of  the  experimental  data  for  those  values  of  femtosecond  TPA  cross  sections 
that  were  measured  under  conditions  similar  to  those  of  our  experiments.  A  numerical  calculation  that  includes  the 
ground  state  absorption  gives  significantly  better  fit  to  the  NLT  data.  These  calculations  have  no  adjustable  parameters. 
A  careful  examination  of  linear  absorption  data  gives  evidence  for  a  two-step  TPA  via  the  Lorentzian  tail  between  the 
S0  and  St  states.  We  conclude  that  the  dominant  contribution  to  the  nonlinear  transmission  in  these  chromophores  in  the 
nanosecond  regime  is  ESA  from  both  singlet  and  triplet  states.  TPA  does  not  contribute  significantly  to  the  transmit¬ 
tance  loss,  but  is  key  to  pumping  the  excited  states.  We  expect  that  this  is  true  for  a  large  variety  of  these  types  of  chro¬ 
mophores.  We  have  also  presented  a  model  of  two-beam  coupling  in  a  nonlinear  absorption  medium  whereby  energy  is 
transferred  from  an  incident  laser  beam  to  an  elastically  backscattered  beam.  In  the  mechanism  proposed,  the  incident 
laser  has  a  nonlinear  time  dependent  phase  and  populates  an  excited  state  of  the  medium  by  two-photon  absorption.  The 
incident  and  backscattered  waves  superpose  to  form  an  interference  pattern,  leading  to  the  formation  of  a  Bragg  grating. 
The  Bragg  grating  consists  of  an  index  modulation  resulting  from  a  modulation  of  the  excited  state  population.  This 
grating  is  out  of  phase  with  the  interference  pattern  which  forms  it  due  to  the  finite  lifetime  of  the  two-photon-populated 
excited  state  and  the  frequency  chirp  of  the  two  waves.  Energy  flows  one-way  from  the  higher  frequency  part  to  the 
lower  frequency  part  of  the  coupled  waves.  Complete  energy  conversion  to  the  backscattered  wave  is  prohibited,  how¬ 
ever,  in  part  because  of  loss  due  to  two-photon  and  excited  state  absorption.  Nevertheless,  the  power  spectrum  of  the 
backscattered  wrave  can  be  nearly  identical  to  that  of  the  incident  wave.  This  mechanism  may  be  at  work  and  explain  the 
NLT  curves  at  lower  intensities  near  energies  of  1 0  pJ. 
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